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For a spatially flat Friedmann model with line element ds 2 = a 2 [da 2 / B(a) — dx 2 — 
dy 2 — dz 2 ], the 00-component of the Einstein field equation reads 8nGToo = 3/a 2 
containing no derivative. For a nonlinear Lagrangian C(R), we obtain a second- 
order differential equation for B instead of the expected fourth-order equation. 
We discuss this equation for the curvature-saturated model proposed by Kleinert 
and Schmidt a. Finally, we argue that asymptotic freedom G~g — •» is fulfilled in 
curvature-saturated gravity. 



1 Introduction 

In the past decades, several extended theories of gravity have been proposed, 
whose effective actions are more general than the Einstein-Hilbert action. 
This approach is motivated by unification schemes which consider gravity 
at the same level as the other interactions of the elementary particles. In 
such theories, we have to define an effective gravitational coupling G e g and 
an effective cosmological constant A e g which give the today observed values 
G e g — > G and A e ff — * A in the weak-energy limit. 

These extended theories introduce new features into gravitational physics 
which general relativity does not posses, in particular higher-order terms in 
curvature invariants or nonminimal couplings between geometry and scalar 
fields. For example, asymptotic freedom can be related to singularity-free 
cosmological models. Its emergence could be the result of thfifact that gravity 
is "induced" by an average effect of the other interactions a t . 

However, the concept of gravitational asymptotic freedom is not com- 
pletely analogous to that in non-abelian gauge theories of strong interaction, 
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since a full quantum theory of gravity does not yet exist. So far, asymptotic 
freedom has been shown for several classes of gravitational Lagrangians which 
are not of physical interest. 

Another interesting feature is that, by extended theories of gravity, cosmo- 
logical singularities can be avoided introducing-a limiting curvature hypothesis 
for values of curvature near the Planck scale This fact is particularly rel- 
evant for string-dilaton gravity since, by duality, it allows to recover large 
classes of cosmological solutions for t — > — oo B. 

In Ref. p, it was argued that the Lagrangian for gravity should re- 
main bounded at large curvature. One example for such a behaviour is the 
curvature-saturated Lagrangian 

£ rs = _J R (i) 

CS 16ttG VI + l 4 R 2 ' 1 ' 

It has been discussed in E ..with ..methods developed for the study of fourth- 
order gravity developed in □ and 13 and the references cited therein. In Eq. (|l|), 
I is a length parameter, and for I — 0, Ccs reduces to the Einstein-Hilbert 
Lagrangian. Earlier models with a non-linear Lagrangian C(R) made only 
polynomial approximation like 

R n 

£ = I^ + E^ ^0 (2) 

k=2 

sometimes accompanied by a logarithmic or a i? m -term with non-integer m. 
In all such cases, one had dC/dR — > ±oo for \R\ — > oo. In contrast to this, 
£cs °f Eq- @) has the behaviour dC/dR — > for R — > ±oo. 

This limiting property can be reformulated in terms of the effective grav- 
itational coupling 



G e ff = 



1D7T-— 

dR 



(3) 



as follows: In the previously extended models Eq. (g), G c g is bounded (in- 
cluding zero) , in contrast to the curvature-saturated model Eq. (Q) where one 
gets |G c ff | — > oo as \R\ ~ > oo; of course, for I — we recover G c ff ^,G. 

Other limiting behaviours have been discussed in □ and □. In 13, a scalar- 
tensor theory with Lagrangian 

C v = F(<p)R+±g*<ps<pj-V(<p) (4) 
has been used where i, j = 0, 1, 2, 3 and the effective gravitational constant is 

Goff = rn^pj (5) 
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instead of our Eq. (||) . In H, G e ft was assumed to depend on the matter density 
in such a way that it vanishes for high density, thus leading to asymptotic 
freedom at high energies. 

In LI, possible finite-size Casimir-effects to the free energy have been cal- 
culated for massive and massless scalar fields, which can produce a quantum- 
effected effective gravitational constant. 

From a geometric point of view, one may argue as follows: For large 
curvatures, the gravitational action / = j Cq^J —gd^x shall not depend on 
any scale; the simplest Lagrangian leading to such a behaviour is 

C c = C m C^ kl , (6) 

the Lagrangian of the conformally invariant Weyl gravity, see i and the ref- 
erences cited there. 

In the present paper, we extend the discussion of Ref. u to more general 
types of curvature-saturated Lagrangians thanEq. (@), and we give more de- 
tails about the set of spatially flat Friedmann models solving the correspond- 
ing field equations. As it will become clear below, the concept of asymptotic 
freedom is meaningful also for curvature-saturated theories and it is necessary 
to define meaningfully the notion of an effective cosmological constant. 

2 Deduction of the Field Equations 

For the action / = J C(R)^J —gd 4 x, it is useful to write the metric of the 
expanding spatially flat Friedmann model as 

T da 2 1 

ds 2 = a 2 —— - dx 2 - dy 2 - dz 2 , (7) 

with a > and B(a) > as shown in 0. In these coordinates, the parameter 
a is called curvature time because the curvature scalar has the simple form 

which is linear in the only unknown function B(a), and it does not contain 
second derivatives. The non-vanishing components of the Christoffel affinity 
are as follows 

r°ao= 1 --^, rl = -e T% = -5 aP , (9) 

a IB a a 

where the dash denotes d/da. The greek indices assume the spatial values 
a, (3 = 1,2,3, and y/—g = a 4 / \f~B. The Ricci tensor has the components 

^(S+Ck* Ro °=^-S' (io) 
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Together with the metric Eq. (Q) , its mixed- variant version can be calculated 
to 

*~(i + 5?)«- 

The most often used direct way to deduce the field equation is to insert these 
expressions into the fourth-order field equation SI/dgij = following from the 
variation of the action /. 

Here we present a shorter and more direct derivation: We first use the 
fact that the field equation also implies 51 /SB = 0, and in a second step, 
verify that in spite of this simplified variation, no spurious solutions appear. 

Let us denote dC/dR by h(R), and d 2 C/dR 2 by k(R). The vacuum field 
equation reads 



° 5B dB V dB ' 

After multiplying this by a 7 to avoid negative a-powers, we obtain 

= -a 7 £(-3B'/a 3 ) + 3a 3 (2B - aB')h{-3B' /a 3 ) + 

+18B(3B' - aB")k(-3B'/a 3 ) . (12) 

It is remarkable that this equation is of second order for one function B only, 
but nevertheless, it is equivalent to the whole fourth-order field equation for 
the metric ([?]). One order reduction follows from Eq. (^), the other from the 
fact that Eq. (^) is a constraint and not the full dynamical equation. 

To exclude the existence of spurious solutions it now suffices to see that 
exactly two free initial conditions can be put: B(ao) and B'(a^). 
Example: Let C = R + 2A, where A is a constant. Then h — 1 and k = 0, 



and Eq. ( |12[) reads 



B(a) = ^. (13) 

The Hubble parameter H is related to B via 

B{a) = H 2 a 4 , (14) 

(see 0). So, with Eq. ( [l3|) we get A = 3H 2 = constant consistent with the 
usual de Sitter space-time calculation. 

3 Solutions of the Field Equations 

Now we consider some more general cases for C(R) and look for the corre- 



sponding solutions of Eq. (12) 
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3.1 Solution for Einstein's theory 

To get a feeling for the coordinates (0) we first look for the 00— component of 
the Einstein field equation 

SirGTij — Rij — -^Rgij ■ 



Using Eqs. (g) and (10), we obtain 

8ttGT 00 = A . (15) 
a z 

The energy density is always non-negative for spatially flat Friedmann mod- 
els in Einstein's theory. Eq. (15) provides us with another justification for 
calling this coordinate curvature-coordinate: a is chosen such that Too does 
not depend at all on the function B(a). 
Raising one index of Eq. @ we find 

8nGp = &nGT° = ^ . (16) 

Together with Eq. (|l4|), this represents a good consistency test: it yields the 
Friedmann equation 

%ixGp = SH 2 , (17) 

i.e., curvature time and synchronous time give identical results. 
For the equation of state p = ap we have 

pa 3{1+a) = constant . (18) 



Together with Eq. (|16[), this yields 

B = B a a 1 - 3a , (19) 
with a positive constant Bq. This is consistent with Eq. (|l3| ) for a = — 1. 

3.2 Solution for C = R m . 

For the Lagrangian C = R m , with constant m (m ^ 0, 1), Eq. ( fl2| ) simplifies 
to 

= 2m(3m - 4)BB' + (m - l)aB 2 

-2m(m- l)aBB" . (20) 

We insert 

B = exp / /3da . (21) 



asympt: submitted to World Scientific on February 7, 2008 5 



The integration constant in Eq. ( |21| ) need not be specified, because Eq. 
implies that with B(a) also CB(a), with any positive constant C, represents 
a solution. 

Eq. @ implies B' = 0B and B" = (/?' + fi 2 )B. Then Eq. pi) becomes 



= 2to(3to - 4)/3 - (to - l)(2m - l)a/3 2 
— 2to(to — l)a/3' . 

Here we insert (3 — 7/a, f3' = j'/a — 7/a 2 and get 

= 2rn(4m - 5)7 - (to - l)(2m - l)a7 2 
— 2to(to — 1)07' . 



By putting x = In a, Eq. (23) can be rewritten in the form 

4m — 5 2to — 1 



dx 



7 



1 



2to 



-1 



(22) 



(23) 



(24) 



The de Sitter space-time is represented by 7 = 4. From Eq. (|24[), it becomes 
clear that it is a solution for m = 2 only, i.e. for L = R 2 . 

All solutions of ( pi[ ) can be given in closed form: For m = 5/4 we get 
with any constant xq- 



7 



3(x - x ) 

and together with ao = e x ° and Eq. ( pl| ) finally 



S(a) = exp 



da 



aln(a/ao) 

5/3 



= Bo In — 

a 



For m = 1/2 we find 7 = ie 6 ^-* ' and 

B(o) = exp[±(a/a ) 6 ] 



(25) 



(26) 



Let us study the stability of these solutions; instead of directly comparing 
with neighboring functions B(a) it is easier to consider the neighborhood of 
the constant solution: 7 = constant in Eq. (]2~I|). Up to the uninteresting 
solution 7 = representing R = we find 



7 = 7o 



4to — 5 2m 
to — 1 2to — 1 ' 
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i.e. 



d-l , , 2m — 1 . 

t x =^-^^t- (27) 

For m > 1/2 this implies stability. For m = 2 this proves again the attractor 
property of the de Sitter space-time. 

In synchronized time, power-law inflation is described by a(t) ~ t n , n > 
1. In our coordinates, such solutions correspond to a constant 7 within the 
interval 2 < 7 < 4. 



3.5 Solutions for Curvature-Saturated Lagrangians 
As another model, consider the high-curvature ansatz 

£ = A+^, (28) 

in which h(R) = -C/R 2 , k(R) = 2C/R 3 , with constants A and C. The 
ansatz (53) is supposed to approximate the regions R — > ±00, and it will be 
matched at small curvature values to a polynomial Lagrangian of type Eq. 
(||). Here we are only interested in the high-curvature regions. The concrete 
values of A and C are not yet fixed, and they may be different for R — > +00 
and R — ► — 00, respectively. 

For A = we can directly apply Eq. ( p7| ) with m = —1, i.e., 70 = 3, and 

7). (») 

This equation can be integrated in closed form, but we only need that 7 = 3 
represents power-law inflation and that this represents an attractor solution. 
In synchronized coordinates, it reads 

ds 2 = dt 2 -t 4 (dx 2 +dy 2 + dz 2 ). (30) 

It gives R — * as t — » 00. Thus, for sufficiently large t, the development leads 
to our Universe today. 

For A 7^ 0, we insert (J^) into (^), and have to solve 

3A(B') 3 = 2a 2 C[5BB' + aB 2 - 2aBB"\ . (31) 
This equation will be integrated further in future work. 
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4 Curvature-Saturated Gravity with Matter 

To include matter, we have to replace the l.h.s. of Eq. (|l2|) by a 7 p. For ideal 
fluid, we can use Eq. ( p"8| ) as an equation of state. As a test of this procedure, 

we can insert C = R/16ttG, and get a 7 p — , which is consistent with 

■ . 8nG 

Eq. ©. 

Therefore, matter with equation of state p = ap can be described by 
p = poa~ 3 ( 1+a ^ and the equations of motion are 

a 7 p = -a 7 £(R) + 3a 3 {2B - aB')h(R) 

+185(35' - aB")k{R) , (32) 
R = -3B'a~ 3 . 

We stress that such equations hold for general £(i?)-Lagrangians and for any 
perfect fluid matter. 



5 Conclusions 

In this paper, we continued the discussion on properties of the curvature- 
saturated cosmological models proposed in Ref. El. Especially, we deduced 
the main equation with matter, Eq. (|32]), for a spatially flat Friedmann model 
in the curvature coordinates (Q). 

We have shown that from the singularity, the Universe expands via power- 
law inflation (which represents a transient attractor) to the actual state. 

In order to find-the Wheeler-de Witt equation for fourth-order gravity 
models (see Refs. O'tSO and the papers cited there), one usually has to 
introduce, in a more or less natural way, more degrees of freedom to reduce 
the equation down to second order. Here we presented a version where the 
Lagrangian can be directly used: C(R)^/—g depends on a, B and B' only; 
of course, differently from the use of the synchronized time coordinate, we 
have now an explicit dependence of the metric coefficients on the time-like 
coordinate a. 

As a final remark, we see that asymptotic freedom G~ s = 16ndC/dR — > 
can easily be incorporated also in curvature-saturated gravity but its meaning 
is different from that in 13. There G e ff is a function of matter density which 
regulates the gravitational coupling, here it is the scalar curvature (i.e. the 
form of the gravitational Lagrangian) which leads towards the saturation and 
then towards asymptotic freedom. Future studies will be devoted to more 
physically motivated effective Lagrangians. 
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